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Abstract 

We report on some implications of the theory of turbulence developed 
by V. Yakhot [V. Yakhot, Phys. Rev. E 57(2) (1998)]. In particular 
we focus on the expression for the scaling exponents We show that 
Yakhot's result contains three well known scaling models as special cases, 
namely K41, K62 and the theory by V. L'vov and I. Procaccia [V. L'vov & 
I. Procaccia, Phys. Rev. E 62(6) (2000)]. The model furthermore yields a 
theoretical justification for the method of extended self-similarity (ESS). 

1 Introduction 

The most challenging open problem of turbulent fluid motion is the phenomenon 
of small scale intermittency, i.e. the unexpected frequent occurence of very large 
fluctuations of the local velocity over small spatial distances. These fluctuations 
on a certain length scale r are usually investigated by means of the longitudinal 
velocity increment u{r), the difference of the velocities at two points in space 
separated by the distance r: u{r) = v{x + r) — v{x) (here, v is the component of 
the velocity fluctuations in direction of the separation vector r). The statistics 
of u{r) is commonly characterized by means of the moments S"{r) = {u{r)"' ), 
the so-called velocity structure functions. 

Even though the first theory of fully developed turbulence dates back to 
1941 [1], a commonly accepted model for the structure functions has yet not 
been found (see [2] for an introduction into the field). The only exception is 
the third order structure function for which an exact relation can be derived 
directly from the Navier-Stokes equation, Kolmogorov's famous four-fifths law. 
For scales r much smaller than the system size L and larger than the dissipation 
length scale of the flow configuration, the third order structure function is a 
linear function of the scale r [2] : 

S\r) = --^er. (1) 
Here, e is the mean rate of energy dissipation within the flow. 
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In the light of the exact result (1), Kolmogorov based his models on the 
hypothesis that the structure functions follow power laws in r 



S'"(r)ocr«", (2) 

where the Cn are the so-called scaling exponents. 

Assuming a cascading process within the turbulent flow which transports 
energy from large to small scales with a constant rate of energy transfer, Kol- 
mogorov [1] predicted a simple linear dependence of the scaling exponents on 
their argument n (furtheron refered to as K41): 

Qn = |. (3) 

However, as noted by Landau [3], there is reason to doubt the assumption of a 
constant rate of energy transfer inbetween scales. Taking into account Landau's 
criticism, Kolmogorov [4] proposed an advanced model for the energy transport 
from large to small scales which, in better agreement with experimental results, 
predicts a quadratic dependence of the scaling exponents on n (this result will 
furtheron be refered to as K62): 

n n{n - 3) 
^" = 3 

For n = 2 it follows from (4) that C2 = | + The parameter 5 thus describes 
a correction to the K41-value of | and will furtheron be denoted as anomaly 
parameter [5]. Its value is not given by the model but has to be determined 
experimentally, the commonly accepted value being 5 = 0.029 ± 0.004 [6]. 

This parameter has gained quite some importance within the theoretical 
framework recently developed by V. L'vov and I. Procaccia [7]. Investigating 
the behaviour of multiscale correlation functions for the case that two or more 
of the scales approach each other, these two authors derive an expression for 
the C,n which is of second order in 8: 

Cn = ^-^^^^.5{l + 2<5(n-2)6}, (5) 

where h is of order —1 with a weak dependence on details of the calculations. 

In the sequel, the theory propsed by V. Yakhot [8] will be discussed with 
respect to these three models. In particular, it is shown that equations (3-5) 
are obtained as the zero, first and second order Taylor expansions of the scaling 
exponents derived by Yakhot. The theory also fixes the parameter b in equation 
(5). Furthermore, it is shown that the theory provides a justification for the 
probably most important experimental technique to measure scaling exponents, 
namely the extended self-similarity (ESS) [11]. 

The paper is organized as follows. Section 2 summarizes some important 
results of Yakhot's theory while section 3 gives the main results of our consid- 
erations. A discussion in section 4 will finally conclude the paper. 
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2 Some results of Yakhot's model 



Following ideas by A.M. Polyakov [9], Yakhot proposed a model of fully de- 
veloped turbulenee in tlw^ limit of very large Reynolds number which is based 
on an analytical treatenient of the Navier Stokes equation [8]. For the case 
of turbulent flows definded by some large scale boundary condition (a typical 
example is free jet turbulence), he derived the following equation for the scale 
dependence of the probabiltiy density function p{u) of the velocity increment: 

Here, A and B are two parameters that are not fixed by the theory and a = 
\J {v'^) is the rms of the velocity fluctuations v (with {v) = 0). 

By multiplication of eq. (6) with and subsequent integration with re- 
spect to u, the equation for the structure functions S"{r) = J u^p{u)du can be 
derived. The resulting equation reads: 

r|:5"(r) = Q^S-{r)+aZn^S-\r), (7) 

= (8) 



B + n' 
n{n — 1) 

B + n 



(9) 



In the limit of small scales r <^ L, the second term on the rhs of equation (7) 
becomes neglible. In this limit, the solutions of the resulting equation are simple 
power laws: 

5"(r) = C„r'^". (10) 

The four-fifths law (1) imposes the condition (^3 = 1 on the scaling exponents 
which, inserted into equation (8), leads to: 



B + 3 



(11) 



The scaling exponents can thus be expressed as a function of only one yet 
unknown parameter: 

Cn = ^^. (12) 

Taking into account both terms on the rhs of eq. (7) , it is furthermore possible 
to determine the parameter B. For the second order structure function the last 
term on the rhs of eq. (7) vanishes on all scales (as (r) = 0) and the general 
solution is a power law in r: 

S\r) = C2r^^ (13) 

At the integral length scale L the velocity fluctuations are uncorrelated. This 
condition fixes the integration constant C2 in eq. (13) and the final solution for 
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the second order structure function is [8]: 

S\r) = 2a^Q'\ (14) 

Inserting this result into the equation for the third order structure function leads 
to: 

= sHr)+2ahs[-) . (15) 
The general solution for this equation is: 

SHr) = C3r + 2.3|(L)^-\ (16) 

The four-fifths law (1) immediately determines the integration constant C3 and 
the final result for the third order structure function is: 

n, , 4 / r\C2+i 

SHr) = -5er + 2a3^(-) 

4 _ „ B + 2 /r\C2+i 



= --er+18c7 



5 (S + 3) 



(z) • (1^) 



The parameter B can be determined from the condition that odd order moments 
vanish at the integral length scale, i.e. S'^lL) = 0. As furthermore eL a^, it 
follows from eq. (17) that 

t « (18) 

which is solved by S « 20. With this value, the scaling exponents (12) are, 
within the errors, indistinguishable from the best experimental values [8]. 



3 A Reformulation 

An interesting result can be derived from Yakhot's theory by replacing the 
parameter B in the above equations by the scaling anomaly 6. In accord with 
equations (4) and (5) we define S as: 

C2 = l+S. (19) 
On the other hand, we have from equation (12) that 
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Using these two relations to replace the parameter B in eq. (12) leads to the 
following expression for the scaling exponents: 

_ n 2 + 2.5 

3 2 + 3(n-2)(5- ^ ^ 

As (5 is a small parameter, it is possible to expand the in a Taylor series 
around 5 = 0. Up to order two in 5 this expansion reads: 



Cn « - |l - -(n - 3)5 +-(n-3)(n- 2)^2 



This is exactly the expression derived by L'vov and Procaccia, equation (5), and 
in agreement with the order of magnitude they find we obtain h = —3/4. 

It is also interesting to reformulate the equation for the third order structure 
function, equation (17), in terms of 6: 

4 3 V /r \C2+i 



.27Sa^(^f^\ (23) 
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In the second line we neglected terms of order 5^ and higher. Setting S^{L) = 
we finally obtain 

4 eL 4 

in full agreement with the experimental value of 5 = 0.029 ± 0.004 [10]. 

Yakhot's theory also allows for an interesting statement on the moments 
T"(r) = (|M(r)|") of the absolute value of u{r). From equation (6) for the pdf 
p{u), the corresponding equation for the pdf p{—u) of the negative increment 
can c;asily be derived. Adding the equations for p{u) and p{—u) and using the 
definitions 

p'^{u) := p{u)+p{-u) 

p^{u) := p{u)-p{-u) (25) 



we obtain: 



^dp^ _ d_ (dp^\ ^ _ ^A(yp+) 



dr du\ dr ) r du 

Ldv? 



+ tI-Aup-). (26) 



As the T^(r) can be expressed by means of p"^: 

T"(r) = j vJ'p+{u)du, (27) 
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the equation for these moments can be obtained by muhiplying equation (26) 
with u" and subsequently integrating with respeet to u. The result is: 

r-^T^r) = CnT"(r) + azn^r-\r), (28) 

where we defined the moments t"(r) as i"(r) = J^^"^ u'^p^ {u)du (these quanti- 
ties can be thought of as generalized odd order moments). 

Equations (7) and (28) for the odd order moments T^"-+^(^r) and S'^"+-'(r) 
differ in the last term on the rhs. While the equation for 5'^""'"^(r) incorporates 
the large even order moment 5^" (r) , the moment t'^" (r) contained in the equa- 
tion for y2ri+i^^-j jg considerably smaller and even tends to zero as r approaches 
the integral length scale. In other words: While this term dominates the be- 
haviour of the 5^"+^(r) for large scales, it has a much weaker influence on the 
moments T^"+^(r). According to this result, odd order moments of absolute 
velocity increments can in a much better approximation and over a wider range 
of scales be described as power laws in r. This finding is in full agreement with 
experimetnal results [11]. Equation (28) furthermore tells us that the scaling 
exponents determined from the moments r^"+^(r) are indeed identical to the 
exponents that determine the scaling of the structure functions S'^"'"'"^(r) in the 
limit of smale scales. Things even improve if, as it is the case for the extended 
self-similarity, the relative scaling of structure functions is considered. Basically, 
this method assumes a linear dc;pendcnce between the logarithms of structure 
functions of various order [11] where the slope is given by the scaling exponents. 
The ESS-property can be formulated as: 

dlnT^jr) ^ Cn dluTPjr) 

dlnr Cp d\nr ' ^ ' 

An equation quite similar to the ESS-property (29) can be derived from Yakhot's 
model. Equation (28) for the T"(r) can be rewritten as 

1 rainT"(r) rr-i(r) 

1 = ^ < TT, CrZr,. — 



C„ i ainr "L T"(r) 

1 \d\nTP{r) r tP-^{r) 



Cp 1 d\nr TP{r) 



(30) 



which finally yields: 



ainT"(r) Cnd\nTP{r) r ( t-^-^r) p-ltP-\r) 

- +(^z„-< - \ . (31) 



ainr Cp d\nr L [ T'^{r) n-1 TP{r 

Apart from the last term on the right hand side, this equation is just the ESS- 
property as expressed by equation (29). It is furthermore found that this last 
term not only inlcudes the small quantities (r)/T" (r) but even the difference 
of two of these small terms. It thus turns out that the extended self-similarity 
minimizes the influence of this term and in fact measures the scaling exponents 
that describe the small-scale structure of turbulence. 
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4 Conclusion and Comments 



From a theoretical point of view the most interesting of our results is the find- 
ing that Yakhot's model unifies three of the most important theories of strong 
turbulence: It not nonly contains K41 and K62 as limiting cases but also the 
scaling exponents predicted by L'vov and Procaccia. Besides, Yakhot's theory 
applies to a much wider range of scales and physical phenomenona: The results 
presented here were obtained by taking two limits, the limit of small scales r 
and the limit of a small anomaly parameter 5 (see [8] for a discussion of the 
extreme cases i? = and B +oo). 

Another important result follows from Yakhot's model for the extended self- 
similarity and the use of moments of absolute velocity increments. The theory 
gives evidence that both methods in fact allow for a precise and reliable deter- 
mination of the scaling exponents that describe the small-scale behaviour of the 
structure functions. 

Finally, we would like to point out that the model equations such as eq. (24) 
predict a dependence of the scaling parameter 5 on the large scale properties 
of the flow. This dependence, though obviously being weak, is in contradiction 
to the commonly assumed universal properties of small scale turbulence. It is, 
however, in agreement with recent experimental findings [12, 13]. 
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